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A pronounced local in-gap zero-energy bound state (ZBS) has been observed by recent scan¬ 
ning tunneling microscopy (STM) experiments on the interstitial Fe impurity (IFI) and its nearest- 
neighboring (nn) sites in FeTeo.sSeo.s superconducting (SC) compound. By introducing a new 
impurity mechanism, the so-called tunneling impurity, and based on the Bogoliubove-de Gennes 
(BDG) equations we investigated the low-lying energy states of the IFI and the underlying Fe- 
plane. We found the peak of ZBS does not shift or split in a magnetic field as long as the tunneling 
parameter between IFI and the Fe-plane is sufficiently small and the Fe-plane is deep in the SC 
state. Our results are in good agreement with the experiments. We also predicted that modulation 
of spin density wave (SDW), or charge density wave (CDW) will suppress the intensity of the ZBS. 

PACS numbers: 74.70.Xa,74.55.+v,74.25.N-,71.20.-b 


Since the discovery of iron-based superconductor, [IJ 
new compounds continue to be found. The 11 type 
iron chalcogenides has attracted much attention due to 
the simplicity of its crystal structure. Angle-resolved 
photoemission spectroscopy, STM, and transport experi¬ 
ments have been performed to investigate the elec¬ 
tronic structure and SC gap. By substituting Se for Te in 
FeTe compound, superconductivity appears and presents 
a variety of phenomena for different mixing ratios of Te 
and Se .6§. The as-grown Fei +y Te x Sei_ x single crystals 
contain a large amount of excess Fe that are randomly 
situated at the interstitial sites in the crystal. Recent 
STM experiments investigated the electronic state near 
the interstitial Fe in FeSeo.sTeo.s superconductors. Q A 
robust ZBS is observed in the tunneling spectrum taken 
at the center of the IFI and its nn sites with the intensity 
decaying quickly. Interestingly, the peak does not split 
or shift in a magnetic field which is drastically inconsis¬ 
tent with the magnetic or non-magnetic impurity effects 
in either d-wave or s-wave superconductors. [8J-fl0j 


The interplay of the IFI with the in-plane Fe and the 
opposite Se (Te) are complex. The valence of excess 
iron does not equal to that of the in-plane iron, the 
effects of the IFI to the physical property are sensitive 
to the stoichiometry. [5, 11-14], Although the magnetic 
moment of the access Fe has finite value, for SC state 
Fei + yTe 0 . 5 Se 0 . 5 , we assume that the IFI is coupled to the 
underlying Fe-plane by a hopping term. Taking the IFI 
as a tunneling impurity {TH} without scattering potential, 
and based on the BDG equations in the presence as well 
as in the absence of the magnetic field, we systematically 
investigated the low-lying energy states of the IFI and 


the in-plane Fe. We found the IFI-induced in-gap ZBS is 
a common feature for iron-based superconductors as long 
as the hopping term between the IFI and the underlying 
Fe-plane is sufficiently small and Fe-plane is deep in the 
SC state. The ZBS is isotropic and localized, the inten¬ 
sity of the induced ZBS decays to zero at a distance of 
three or four lattice constants away from the IFI. 


In the optimally doped regime, and when the external 
magnetic field is not very large, beside the appearance 
of the in-gap resonance peaks induced by the magnetic 
field [16|, 17], the IFI-induced no-splitting no-shifting ZBS 
shows up at its nn sites no matter where we put the 
IFI. This phenomenon is contrary to the Zeeman effect 
of single electron since the IFI-induced ZBS results from 
Cooper pairs. Intensity of the ZBS at nn sites has the 
maximum value as the IFI is located above the vortex 
core center. Considerably large magnetic field may in¬ 
duce charge oscillation, and breaks the process of reform¬ 
ing Cooper pairs when electron tunnels back to the Fe- 
plane from the IFI. Therefore larger magnetic field will 
suppress the ZBS at the nn sites of the IFI. We believe 
that the tunneling impurity is a new scenario to under¬ 
stand the experimental finding of IFI effect. 


Since the electronic structure and Fermi-surface topol¬ 
ogy of the 11 system is very similar to those of the other 
iron-based superconductors, Q we use a minimal two- 
orbital four-band tight binding model [Hj to investigate 
Fei+ y Teo. 5 Seo. 5 . Different from the co-planar Cu02, su¬ 
perconducting Fe-plane is sandwiched between two adja¬ 
cent Se(Te) layers, the next-nearest-neighbor (nnn) hop¬ 
pings are mediated by up and down Se (Te) and are not 
equal to each other. This asymmetry should be correct 









2 



FIG. 1: (color online)Schematic plot of the tight-binding 
model. The right picture shows the relative position of IFI to 
the Fe-plane. 


when one investigates the surface properties as in STM 
experiments, because the bonds between up Te (Se) ions 
and Fe ions are broken after cleavage. Theoretical re¬ 
sults 16, 17, 3-2^ based on this model are qualitatively 
consistent with experiments. [24-26] The projected up 
and down Se (Te) atoms stay at the centers of the Fe pla- 
quette alternatively, and the IFI is located at the opposite 
side of the down Se (Te) above the Fe-plane[see FigQ]. 
The total Hamiltonian is expressed as H = Ho + H\fi, 
Hq is the mean-field phenomenological model Hamilto¬ 
nian without excess Fe, reads 


Ho — y (tub i' c ij.a c ji'cr + h.c.) ji y 

i/xjz/<7 i ii(T 

+ U 'y \ni)nj vcr + + h.c.) , ( 1 ) 

i/xj vcr 

where n^a is the electron concentration at site i, orbital /i 
for spin <r, y is the chemical potential, determined by the 
average electron filling on Fe site (n)i = 2 + x, x denotes 
electron doping concentration. Hq can be separated to 
three parts. The first line is tight bindin g p art, with the 
hopping integral fi _4 = 1,0.4,—2,0.04 [18| as depicted 
in Fig [IJ The first term in the second line is the mean 
field expression of Coulomb and Hund’s interaction Ho.mt 
in the SC state without SDW or CDW. It will be more 
complicated when SDW or CDW exists 

Ho^int \ E T (U SJh') {P'iu' 

i ,v'^v,cr 

+ {U -2J H ){n- 1V 'a)n- lva . (2) 

Here we take U = 3.4, Jh = 1.3 and U' = 3t/ ~ 4 5l7g . [l9l — 
[23j The last term in the second line is the phenomeno¬ 
logical superconducting pairing part. Although there 
are controversies about the symmetry of the SC pairing, 
s_l_-wave 3, mil is suggested for Fei +y Tei_ x Se x . In 
real space it is intraorbital pairing = w — 

where j is the nnn of i site. The self-consistent 
mean field are Ai^ and (n i M ). 

The role of the IFI is to provide two tunneling channels 
corresponding to the two orbitals of Fe. Electron tunnel¬ 
ing onto IFI belongs to a Cooper pair since the Fe-plane 


is in the SC state, when it tunnels back to the Fe-plane a 
new Cooper pair reforms. Because the effective magnetic 
moment of the cooper pair is zero, there is no magnetic 
interaction between the IFI and the Fe-plane. The effec¬ 
tive coupling of the IFI to the Fe-plane reads 

Hifi = —ti fi y C^Ci^cr , (3) 

(i/tcr) 


where c 1 is the creation operator of the IFI, () means 
the summation up to its four nn sites in the Fe-plane. 
Tunneling magnitude Hfi is a tuning parameter in our 
calculations. 

Magnetic field contains a dynamical term for electrons, 
in the mixed state, the effect of the magnetic field is in¬ 
cluded through the Peierls phase factor. For a perpen¬ 
dicular magnetic field, the hopping integral in Eq. © and 
© should be changed to t' which can be expressed as 


exp [i <pij], where ^ f- A(r) • dr, and 

$0 = he/ 2e is the superconducting flux quantum. We 
assume the applied magnetic field B to be uniform and 
the vector potential is A = (—By, 0, 0) in the Landau 
gauge. Periodic boundary condition should be consistent 
with the requirement that the total phase factor along 


each small plaquette given by ipij = . 

The presence of Peierls phase makes the usual transla¬ 
tion operator change to magnetic translation operator. 
To ensure them commutable with each other and the 
Hamiltonian, each magnetic unit cell has to contain 24>o 
flux 29, 3fJ. In this case the supercell technique can be 
used to calculate the local density of states (LDOS) in 
the presence of B. We self-consistently solve the BDG 
equations in real space, linear dimension of a unit cell is 
N x x N y = 32 x 32 for B = 0, while for finite B we adopt 
N x = 2N y . The LDOS is calculated according to 


PiM = y [KVkP'K-En.k - w) + \v^\ 2 S(E nM + w)], 

n/ik 


where the delta function <5(a;) is r/7r(cc 2 + T 2 ), with the 
quasiparticle damping T = 0.005. and are the 
eigenstate component of the energy E n ,k at site i orbital 
/r for wave vector k, and they correspond to the particle¬ 
like component and hole-like component respectively. 
The number of unit cell is taken as M x x M y = 20 x 20 
for B = 0 and 10 x 20 when external B is applied, with 
k Q = — —|- , a = x,y. Throughout the paper, 

temperature is zero K, the energy and length are mea¬ 
sured in units of ti and the nearest Fe-Fe distance a, 
respectively. 

Although the experiments was performed on half-filled 
SC compounds Q, we found ZBS is a common feature 
for SC state as long as Hfi is small. Figj2] shows LDOS 
p for different fiFi with x = 0.0 and x = 0.12 for the 
Hamiltonian of Eq[T| For x = 0.0 and Hfi less than 0.04, 
LDOS at the IFI has a single sharp ZBS peak which can 
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FIG. 2: (color online) LDOS at IFI and its nn sites as 
well as at the nearby nn + 2 site for different tipi = 
O.Olti, 0.02ti, 0.03ti. Panel (a)-(c) are for x = 0.0 while panel 
(d)-(f) for x = 0.12. In (d) and (e) the dotted volet line are 
for tiFi = 0.05. 


be seen in Figl^a), and the peak height is decreased with 
the increasing of fiFi • 

At its four nn sites p is isotropic. Compared to the 
pure sample, a pronounced in-gap ZBS appears and the 
height of the peak is enhanced with tipi increasing from 
0.01 to 0.03 as seen in Figl^b). At the sites two lattice 
constants away from the nn sites, the intensity of the 
in-gap ZBS decays almost to zero as shown in Figl^fc). 
For different doping cases, we expect the ZBS still exist 
whenever Fe-plane is in the SC state. Fig[2])d)-(f) plot 
p at the IFI and its vicinity for x = 0.12. The curves 
are very similar to those of the x = 0.0 cases with the 
corresponding peak of the ZBS a little higher. 

Small and large Iifi give significantly distinct results. 
For relatively larger tipi, the results deviate from the ex¬ 
perimental observation. When x = 0.12, tipi = 0.05, p at 
the IFI splits into two asymmetrical peaks shown by the 
violet dotted line in Fig[2])d). Due to the proximity effect 
the LDOS at the nn sites of the IFI also shows double 
peaks with heavy weight on the positive energy and can 
be seen in Fig^e). Therefore, sufficiently small tipi is an 
important condition for the appearance of ZBS. As tem¬ 
perature increases, the phonon mediated layer tunneling 
Iifi is also increased, thus ZBS will not appear for higher 
temperature cases. 

We have addressed the cases when Fe-plane does not 
have magnetic order defined as mi = ( n 4 — 714 )/2. In 
fact, as magnetic order coexists with the SC order which 
will happen in underdoped cases, the LDOS at the IFI 
and its nn sites will still have the in-gap ZBS. In the fol¬ 
lowing we will discuss the more complicated cases with 
Ho'int taking the form of Eq{2] The formation of the ZBS 



FIG. 3: (color online) For Iifi = 0.02, x = 0.06, (a)LDOS 
at nn sites without B. (b)LDOS at nn site for B ~ 22T with 
IFI projected at vortex core center (12.5,12.5). (c)Imagine 
of modulated \rm\ for B ~ 22T with IFI above vortex core 
center. 


at the nn sites is similar to the Andreev reflection and we 
expect the homogeneous magnetic order itself does not 
suppress ZBS as long as the sample is in the SC state. 
We can see from Fig|3][a) that the ZBS still appears at 
nn sites of the IFI for underdoped case x = 0.06 with 
Iifi = 0.02. While an external B could drive SDW mod¬ 
ulation [see Figj3][c)] and CDW in this case, the competi¬ 
tion between the oscillated SDW and SC order ruins the 
proximity effect. Therefore as SDW exists, an external 
magnetic field will suppress the ZBS as well as other in¬ 
gap states which is depicted in Fig[3](b). In Figj3][b), the 
IFI is located above the vortex core center with the ZBS 
almost vanished; as the above IFI moves away from the 
center, the ZBS does not appear. 

In the optimally doped regime, we investigated pertur¬ 
bation of IFI on the mixed state, p is plotted in FigQ] 
at IFI as well as at its nn sites for different B with fixed 
Iifi = 0.02, x = 0.12. A unit cell of linear dimension 
64 x 32 corresponds to a relatively small magnetic field 
B ss 13T and 40 x 20 corresponds to relatively large field 
B « 32T. In all cases of small tipi, p on the IFI has a 
sharp ZBS regardless of the strength of B and the posi¬ 
tion of the IFI. The height of the peak becomes lower as 
the IFI moves close to the vortex core center which can 
be seen in Fig|4](a). 

Now we look at the relatively small magnetic field 
B « 13T. In the absence of the IFI, two in-gap resonance 
peaks induced by B are located at negative energies in 
the vortex core. At sites away from the vortex, resonance 
peaks are suppressed and move to the gap edge, finally 
evolve into its bulk feature [see Fig|4](b)]. When intro¬ 
ducing IFI into the system, the above characteristics do 
not change except for the additional appearance of ZBS 
induced by the IFI which is depicted in Figj4][c). We can 
see that the in-gap ZBS shows up at nn sites no matter 
where we put the IFI. Compared with the cases with¬ 
out B [Fig[2], the height of the ZBS at nn sites of the 
IFI is obviously enhanced, especially when the IFI is lo¬ 
cated above the vortex core center and has its maximum 
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FIG. 4: (color online) (a)LDOS at different IFI for 

B ~ 13T.(b)Without IFI, LDOS at vortex core center and 
at sites away from vortex for B ~ 13T.(c)With the IFI lo¬ 
cated at different place, LDOS at the corresponding nn site 
for B ~ 13T. (d)Similar to (c) but for B ~ 32T. From up 
to down the curves correspond to the projected IFI at vor¬ 
tex core center, and L lattice constants away from the center 
along horizontal axis. The gray dashed line indicate the po¬ 
sition of zero energy. (e)Image of rii for B ~ 13T with the 
projected IFI at (16.5,16.5). (f)Image of m for B ~ 32T with 
the projected IFI at (29.5, 15.5). The curves of (b)(c)and (d) 
are displaced vertically by 0.5 unit for clarify. 


value. As the IFI moves away from the vortex center, 
the peak on nn sites is decreased. When the projected 
point of IFI on Fe-plane is 16 lattice constants away from 
vortex core center along x-axis, the height of the ZBS of 
nn sites evolves back to the value of B = 0 case. We 
also notice that a lower height of ZBS on the IFI [see 
FigH[a)] corresponds to a higher peak of ZBS on its nn 
sites [see FigQJc)]. The more of quasiparticle tunneling 
onto nn site, the less part is remained on the IFI. Results 
of B « 22T (48 x 24 unit cell) are very similar to those 
of B « 13T and we do not show them here. 

Without the IFI, m is enhanced and the pairing pa¬ 
rameter A i = is almost vanished at the vortex 

core center, with A,; and nt being symmetric with respect 
to the two vortex cores. Introducing of the IFI breaks the 
symmetry. For smaller B « 13T and projected IFI at vor¬ 
tex core center (16.5,16.5), although m on the two vortex 
core are different, it increases to bulk value at the scale of 
coherence length. Out of vortex core, rii is almost homo¬ 
geneous which can be seen in the FigQJe). While larger B 
leads to oscillation of rii , Fig|4]T) displays the inhomoge¬ 
neous rii for B « 32T with IFI located above vortex core 
center (29.5,15.5). The difference of rii on the two vor¬ 
tex cores is 0.01, much larger than B = 0 cases in which 
Arii induced by IFI is at the order of 10~ 3 . Obviously, 
magnetic field significantly enhanced the impact of IFI to 
the Fe-plane. The oscillated CDW ruins the Cooper pair 
tunneling and reforming, hence suppresses the ZBS on nn 


sites. For large B « 32T, Figj4][d) shows that when IFI 
locates above vertex core center, ZBS appears at its nn 
sites with the peak much lower than that of small B case. 
As the above IFI moves away from vortex core center, it 
still have impact on nn sites, however the induced state 
is not ZBS and almost invisible in Fig@[d). 

In summary, motivated by the recent STM experi¬ 
ments |3] on iron-based superconductor FeTeo.sSeo.s, we 
investigated the impact of IFI on the s±-wave iron-based 
superconductors. Taking IFI as a tunneling impurity, [15] 
we calculated LDOS on the IFI as well as on Fe-plane in 
the absence and in the presence of magnetic field. In all 
cases with and without external magnetic field, LDOS at 
the IFI has a single sharp ZBS as long as hopping pa¬ 
rameter tipi is sufficiently small and Fe-plane is in SC 
state. The ZBS always appears on its nn sites due to the 
proximity effect in the absence of B. When a magnetic 
held is applied, our results are in good agreement with 
the experiments. When Fe-plane is deep in SC state, the 
IFI-induced ZBS shows up at nn sites no matter where 
we put the IFI. It does not split or shift and the intensity 
is enhanced. The effect of the IFI is isotropic and local¬ 
ized in all cases. However when Fe-plane has oscillated 
SDW or CDW, the intensity of ZBS on nn sites will be 
suppressed. In underdoped cases, a magnetic held can 
lead to modulation of SDW. And when the applied mag¬ 
netic held is considerably large, it may induce CDW in 
the IFI system. 

Tunneling impurity only provide tunneling channels 
in the system. Due to the superconducting state of 
Fe-plane, the effective tunneling is Cooper pair tunnel¬ 
ing with zero magnetic moment thus no Zeeman split¬ 
ting. Although experiments and density functional the¬ 
ory study verihed the access Fe is spin polarized (l3[. we 
do not consider exchange interaction between excess Fe 
and superconducting Fe-plane. The origin of IFI-induced 
ZBS on nn sites is related to Andreev reflection. As far 
as we know, a convinced explanation for the STM ex¬ 
periments is still lacking. Although it may associated 
with non-trivial topological or other factors, our simple 
tunneling IFI gives a natural explanation for the exper¬ 
iments and provides a new interpretation for interstitial 
impurity. 
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